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COMPUTATION  OK  VELOCITY  AND  PRESSURE  VARIATION  ACROSS  AXI SYMMETRIC  TH I CK | TURBULENT  STERN  FLOWS 


By 

Thoaas  T.  Huang  and  Ming* Shun  Chang 
David  W,  Taylor  Naval  Ship  Research  and  Development  Center 
Bethetda,  Maryland  20084 


Abstract 

When  the  curvature  of  a  ships  surface  Is 
large  and  the  stern  boundary  layer  is  thick  the 
cross-stream  inviscld  velocity  and  pressure 
variation  becomes  important  in  stern  boundary 
layer  computations.  The  second  momentum 
equation  is  of  the  form 


where  *  is  the  curvature  of  the  stern  flow,  which 
is  different  from  the  surface  curvature  of  the  body, 
u  is  the  velocity  along  a  streamline,  p  is  the 
pressure,  and  p  is  the  mass  density  of  the  fluid. 

Two  viscous- invisc  id  interaction  computation  proce¬ 
dures  are  presented.  One  uses  a  marching  technique 
in  a  natural  streamline  coordinate  system  together 
with  the  k-f.  turbulence  model  to  solve  the  axi- 
symmetrlc  Reynolds-averaged  parabolized 
Navler-Stokes  equations  for  the  streamline  velo¬ 
city  u,  flow  angle  a,  curvature  *c,  and  hence  the 
pressure  p.  The  numerical  procedure  starts  at  a 
stat Lon  on  the  body  where  the  boundary  layer  is 
thin  and  marches  downstream  Into  the  wake.  The 
boundary  conditions  for  the  values  of  u  and  p 
over  the  Inlet  plane  and  along  a  cylindrical 
stream  surface  outside  the  boundary  layer/wake 
are  set  by  the  appropriate  values  obtained  from 
the  other  simple  viscous- inviscld  interaction 
computations  using  a  modified  thin  boundary  layer 
method  and  potential  flow  calculations  about  an 
equivalent  displacement  body.  Comparisons  are 
made  between  the  numerical  results  and  the 
experimental  data  for  four  different  sterns. 
Comparisons  of  measured  axial  and  radial  velocity 
and  pressure  distributions  and  those  computed  by 
the  simple  Interaction  approximations  and  partially 
parabolized  techniques  have  been  made.  The  simple 
and  efficient  viscous- Inviscld  procedure  for  the 
computation  of  velocity  and  pressure  variations 
across  thick  turbulent  stern  flows  has  been  shown 
to  be  accurate  enough  to  use  as  a  prel iminary 
des  i gn  tool . 

K  Int rodoct ion 

Manv  propellers  and  appendages  are  located 
Inside  of  ship  stern  boundary  lavers.  Therefore, 
it  is  essential  for  naval  designers  to  obtain  a  fun¬ 
damental  understanding  and  an  accurate  predictions 
of  tins  special  class  of  external  thick  turbulent 
stern  flows.  A  series  of  experiments  has  been 
conducted  at  David  W.  Tavlor  Naval  Ship  RM>  ('enter 
to  determine  the  unique  turbulence  structure  and 
v  iscous- tnv  isc  Id  interaction  of  thick  ax  isymmet r ic 
1,2,  i, 4}  and  simple  three-dimensional  1*1, b, 7  f  stern 
flows.  The  l.ighthill  1 8  1  displacement  body  concept 
has  been  proven  experimentally  to  he  an  accurate 
approach  for  comput  ing  viscous- invisc id  stern  flow 
interaction.  The  measured  static  pressure  distri¬ 
butions  on  t  lie  bodv  and  across  the  entire  boundary 
lavers  were  predicted  bv  t he  d  I splacement  bodv 
method  to  an  accuracy  within  one  percent  of  dynamic 
pressure. 


Neither  the  measured  values  of  eddy  visco¬ 
sity  nor  the  mixing  length  were  found  to  be  pro¬ 
portional  to  the  local  displacement  thickness  or 
the  local  boundary- layer  thickness  of  the  thick 
axlsymmetric  boundary  layer.  The  measured  mixing 
length  of  the  thick  stern  boundary  layer  was 
found  to  be  proportional  to  the  square  root  of 
the  local  cross-sect  ion  area  of  the  turbulence 
region  (9].  Thi9  simple  similarity  hypothesis 
for  the  mixing  length  and  the  displacement  body 
concept  has  been  incorporated  into  the  Douglas 
C-S  differential  boundary-layer  method  I  ID1  bv 
Wang  and  Huang  {111.  The  method  predicts  satis¬ 
factorily  the  measured  mean  velocity  distributions 
for  the  boundary  layer  flows  of  five  sets  widely 
different  axlsymmetric  bodv  shapes  and  has  neon 
used  as  a  reliable  design  tool. 

Nakayama,  Patel,  and  (.andwcher  (12,  M', 
and  Dyne  {14}  do  not  use  the  displacement  bodv 
method  to  solve  the  interaction  problem.  In  all 
of  the  methods,  the  flow  field  is  divided  into 
an  inner  viscous  region  composed  of  all  or  part 
of  the  flow  in  the  body/wake  domain  and  an  outer 
potential  flow  region.  Differences  arise  in  t In¬ 
equations  used  to  solve  the  viscous  flow  and  the 
manner  of  defining  the  inner  and  outer  regions. 

Dyne  uses  the  integral  approach  of  Head  115!  to 
calculate  the  boundary  layer  flow  over  the  forebody. 
Over  the  stern  of  the  body  and  in  the  near  wake, 
he  uses  a  differential  approach  which  approximately 
accounts  for  the  curvature  and  normal  pressure 
variation  effects.  An  important  feature  of  his 
approach  is  that  the  boundary  layer  equations  are 
solved  alcng  streamlines,  which  leads  to  a  simpl I- 
ficatlon  in  the  form  of  the  equation.  Also,  no 
distinction  needs  to  be  made  between  flow  over 
the  body  and  in  the  wake.  In  the  approach  bv 
Nakayama,  Patel,  and  I^indweber,  integral  relations 
involving  conservation  of  momentum  and  continuity 
are  used  to  relate  the  momentum  area  and  stream 
function  at  the  hodv-wake  junction.  Patel  and  Lee 
(lb)  present  some  results  for  a  differential 
approach  which  includes  all  of  the  curvature  effects. 
A  modified  d  i  spl acement -hod v  method  based  on  u 
simple  pressure  mapping  has  been  applied  to  axi- 
svmmetric  bodies  hv  Hoffman  M7‘.  The  part  iallv 
par.'ibol  ic  flow  assumpt  ions  have  been  used  to 
solve  axlsymmetric  flow  problems  bv  Muraoka  ,18- 
and  Uhen  and  Patel  {19}  and  for  three-dimensional 
flow  past  surface  ships  hy  (lien  and  Patel  il9*  and 
Muraoka  (2()}.  An  efficient  st  reaml  ine- it  erat  ion 
method  with  a  two-equation  k->  turbulence  model  has 
been  developed  bv  Zhou  (21'  to  compute  turbulent 
Mow  around  ax  isvmmet  r  ic  sterns. 

A  numerical  method  using  a  partially  purabo- 
1 ic  marching  technique  in  a  streamline  coordinate 


i 


system  and  the  k-f  turbulence  model  has  been 
developed  by  Hogan  122}  to  compute  the  turbulent 
flow  at  the  stern  and  the  wake  of  bodies  of  revo¬ 
lution.  This  method  uses  the  modified  Douglas 
C-S  method  (3,111  to  provide  the  initial  conditions 
at  the  upstream  plane  and  the  boundary  conditions 
at  a  large  distances  outside  of  the  boundary 
layer . 

In  this  paper,  the  simple  modification  of 
the  Douglas  C-S  computation  method  (10)  will  be 
updated  and  summarized.  Hogan* s  (221  partially 
parabolic  marching  technique  and  the  k-f.  tur¬ 
bulence  model  will  also  be  reviewed  and  improved. 

The  cross-stream  pressure  distributions  and  mean 
velocity  distributions  computed  by  these  two  methods 
will  be  compared  with  the  experimental  data. 

11.  S  imple  Viscous-lnvlsc id  Interaction  Method 

for  Ax isymmetr ic  Stern  Flows 

This  method  is  an  updated  version  of  the 
method  described  in  References  3  and  11,  The  Douglas 
C-S  method  110}  consists  of  using  Keller's  box 
scheme  to  solve  the  following  set  of  partial  dif¬ 
ferential  equations  expressing  conservation  of 
momentum  and  continuity. 


-  u"v"  •  v  ill  •  i  VinY*  for  outer  region,  yr  <  y 

°^y  0  \Jy) 

\>0  -  0.0148  Ytr/“  (Ue-u)dy  -  0.168  UeS*Ytr  O) 

Ytr  -  [l  +  5.5  (y/£)6]-1 

-  Klebanoff's  intermlttency  factor 


0.169  ^y(ro+0-6«2-ro2 
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the  mixing  length  in  the  outer  region  of  the 
thick  boundary  layer,  6^0.23rQ 


°-‘t.  ‘ik) 


(9) 


the  mixing  length  In  the  Inner  region  of 
thin  and  thick  boundary  layer. 


A  =  26v 


,  damping  length 


u 


flu 

3  s 


v 


flu 

fly 


I  in  _  1  iL 

p  As  r  fly 


r(v 


flu 

ly 


u  "v  ' 


>] 


(1) 


4  (ru)  +  h  (rv)  =  0  (2) 

where  u  and  v  are  the  mean  velocity  components 
in  the  s  and  y  directions,  respectively 


s,y  are  the  coordinates  parallel  and  normal 
to  the  body  meridian,  respectively, 


=  r6  /\  dy  »  displacement  thickness 

P  *'°  \  UeJ  (planar  definition) 

6  *  S995  ,  boundary  layer  thick¬ 

ness  where  u/Ufi  = 

0.995 

Ue  is  the  inviscid  (edge)  velocity  used  in  the 
thin  boundary  layer  calculations 

T w  is  the  wall  shear  stress 

yc  is  the  value  of  y  at  which  Vj  -  vQ 


i-  is  the  fluid  density 
p  Is  the  pressure  on  the  body 
r  *  rQ  (s,n)  +  y  cos  a 
rc  is  the  body  radius 
-1  /  drQ\ 

n°  - tan  (dW 

x  is  the  axial  distance  measured  from  the 

nose 

v  is  the  kinematic  viscosity  of  the  fluid 

u\v'  are  velocity  fluctuations  In  the  s  and  y 
directions,  respectively 

TTT  Is  the  Reynolds  stress 


The  flow  in  the  wake  is  modeled  by  the 
following  differential  equation  for  momentum, 
which  is  simply  the  boundary-layer  equation  with 
skin  friction  neglected 


4^  +  (h+2)  ^-4^  -  0 

dx  '  Up  dx 


(5) 


f‘(' -k)k  M 

,  momentum  area 

/=('  ‘  “~,)rdr 

,  displacement 
area 

A/S! 

,  ax isymmet r ic 
shape  factor 

Granville  i 2  3  )  proposes  the  following  equation 
relating  h  to  U0 


The  above  equations  are  the  standard  thin  boun¬ 
dary  layer  equations  with  the  addition  of  the 
transverse  curvature  effect,  where  r  replaces  the 
bod v  radius  r0.  Effects  due  to  longitudinal  cur¬ 
vature  k  and  pressure  variation  across  the 
boundary  layer  are  neglected. 

The  Reynolds  stress  u  v  is  modeled  by 

fly  =  "  for  inner  region, 

Yv  *\r0JUJ  0  -  y  ?  vc 
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where  the  subscript  t  denotes  conditions  at  the 
tail  and  Q  Is  a  variable  coefficient,  for  which 
Granville  recommends  an  average  value  of  7. 

A  measure  of  the  viscous  mass-flux  deficit 
in  the  thick  ax isymmet r ic  boundary  layer  Is 
def  ined  by 

* 

*<ru)  rdr (7> 


/•ro'f"t,r 


Ux(r0)A  */“ 


U  V 


where  ro  is  the  local  body  radius,  6r  the  axl- 
symmetric  displacement  thickness,  <5  the  boundary 
layer  thickness  at  which  ux(r)/Ux(r)  =  0.995,  ux 
the  axial  viscous  velocity,  Ux  the  axial  invise  id 
velocity,  r  *  rQ+ycos'i0,  and  ux  -  ucosao-vslnu0. 
Tims  the  ax isymmet r ic  displacement  thickness 
def  ined  in  equat ion  (7)  becomes 


the  final  solution  is  taken  to  be  the  average  of 
the  values  given  by  the  last  two  iterations. 

The  method  of  using  a  blowing  velocity  distri¬ 
bution  is  used  for  the  final  potential  flow 
computation.  The  lnviscid  velocity  components 
and  pressure  coefficients  across  the  entire 
thick  boundary  layer  are  then  computed. 


if  the  variation  of  lnviscid  velocity  Ux  (r) 

.»*  ross  the  boundary  layer  is  assumed  negligible, 
following  light  hill's  '8:  derivation,  the  boun¬ 
dary  layer  displacement  effect  can  be  represented 
bv  a  source  distribution  on  the  surface  of  the 
hodv  with  its  source  strength  m  of 

2  "  ro  m  '  dT  [2  1,AUx(x’ro)] 

"r  m  ■  Tr-k\[<r°  +  >>'  ■  r°]V**'o>)  (9) 

fhe  <  •  s  |»i  ival  ent  Mowing  velocity  on  the  body  is 
t  hen 

’  >r„  :Uf<r"  +  V2  -  r»]U*(*-ro)|  (l0) 

The  boundary  Inver  equations  (1)  and  (2)  with 
the  modified  mixing  length  and  eddv  viscosity 
:>>r  the  outer  region  of  a  thick  boundary  Inver  {9* 

,  iven  in  equat  ion  (  i)  are  used  to  solve  the  mean 

velocity  component s  u  and  v.  The  axial  and 

radial  velocity  components  are  resolved  bv 
u>; 'uco>,  t-vs  in  t  and  vr  -  usin<  +  vrosi,  re spec - 
i  ivelv.  The  displacement  hodv,  the  corresponding 
.»  iir«e,  and  the  blowing  wlocitv  are  computet) 
a-  <  ord  ing  to  Equations  (H)  ,  (9),  and  (10). 

Tin-  momentum  area  of  the  far  wake  0  can  be 
■  hr  .Trained  hy  equating  the  net  rate  of  momentum 
loss  oj  the  flow  to  the  tot.il  drag  on  the  body. 

!!<•  displacement  area  A  can  In*  determined  bv 
t-qnat  ions  (V)  and  ( M  in  terms  of  the  conditions 

'  'f  and  ht  given  it  the  tail  and  ’0  (h-1  at  x 

if  t  l:«-  tar  wake.  In  tin*  near  wake  region,  where 
it  her  the  houiulaf’  laver  equations  (1)  and  (2) 
i  •  *  :mple  wiki  model  equations  (5)  and  (/>) 

>!■  •  •  .  r  tie,  i  |  itth-degre>  polynomial  is  used 

I-  f  flie  up*sf  ream  md  downstream  <lispl;ic<  - 
‘  '  i  ■  *  i  •  . ,  :  ai  1 1  i  .  t  he  mat  c  h  ing  p«>  int  ■>  X/  1 

■  •  •  f  ■  j  .  >  i  j :•■ }  .  . 1  >5  .  H< »wever  ,  t  in • 

?  »•  hing  point  must  he  moved  upstream 

*  *'■  apar.it  ion  point  whenever  f’ow  separation 

i  '-*  'aii  in.-  ot  the  il  i  spl  ;u  em.-nt  surface 

*  •  '  •  'f  wake  i  ,  M  <•  shortcoming  of  this  simple 
-  it”-  i  .  i  il  mt  i  a<  t  ion  procedure,  which 
>  «  '*  am  m«  !  in  the  light  «•!  a  more  accurate 

o'  I  1  .  .  ♦  m  the  next  sec  t  ion . 


One  obvious  defect  of  the  thin  boundarv-1 aver 
Equations  (l)  and  (2)  is  the  assumption  of  using 
a  constant  pressure  across  the  boundary  layer  for 
large  values  of  v,  where  l  lie  velocity  component 
parallel  to  the  body  tangent  is  equal  to 
U0cos«t0.  However,  the  computed  value  of  u 
approaches  the  lnviscid  velocity  on  the  hodv  l't. 
Instead  of  U  ros**o.  On  the  basis  of  a  large 
collection  of  stern  boundary-layer  data  f,ll  , 
the  lnviscid  influences  on  the  c  imputed  tangen¬ 
tial  velocity  u  is  adjusted  bv 

"m  =  u(y)  lip.osf  m-.p  ,  f-  nn 


where  is  the  improved  tangential  velocity  for  a 
thick  boundary  layer,  u/U^  •  f  is  the  nondimen- 
t  ional  tangential  velocity  predicted  bv  the  thin 
boundary-layer  equations  (1)  and  (2),  and  llQ  is 
the  f roe-stream  velocity,  ij  -  tan  "l  ( x) /^x(  x)  I , 

Tp  =  0,r"+,V^1/2»  aru*  *’p »  I'r .  \  are  the 
radial  and  axial  lnviscid  velocities  calculated  by 
the  final  iteration  of  the  potent  la 1 -flow  computa¬ 
tion  using  the  blowing  velocity  distribution.  Equa 
t ion  (11)  shows  that  near  the  hodv  where  y  t 

t0.  Up  l't, ,  and  ct^  it.  At  the  edge  of  the  boundary 
laver,  where  u/l’c  =  f  =  1  and 

Um  _  Uscos(x0-Q] ) 

Uo 

the  subscript  >-  denotes  quantities  at  the 
edge  of  the  boundary  layer.  As  y  — ,  Ip  *  U0, 
j  -*0,  and  um/U0  *  cosiQ.  Thus,  the  modified 
tangential  velocity  has  the  proper  asymptotic 
value  Uo  cosu  far  from  the  body.  This  simple 
adjustment  is  an  improvement  over  thin  boundary 
layer  theory  where  u/U0  *Up/U0  as  y 


The  corrected  normal  velocity  can  be 
obtained  from  the  continuity  Equat  ion  The 

adjustment  of  u  using  formula  (II)  is  made  in  the 
normal  y-direotion  and  Ins  little  effect  on  tin- 


variation  of  u  in  the  s-direction,  i.c.,  '(rum)/-s 

•(ru)/n.s.  It  follows  f  ri  :n  Kquat  ion  (.')  and  tin 
boundary  condition  that  u  =  u,n-v*-0  at  v-0.  Ihu;,  one 


finds  that  no  adjustment  of  tin*  normal  velocity 


required,  where  v  is  the  normal  velocity  calculated 


hv  substituting  u  into  Equat  ion  (.’)  and  is  "inch 
smaller  than  the  value  ot  u. 


The  axial  and  radial  velocities  uv  and  v 
... {justed  for  the  thick  stern  houndar”  laver  arc 
given  bv 


it*  i  at  ;-n  process  consisting  of  cal<  u-  u  (|  v 

’  it  in,'  pre  -.urc  di--.tr  ibut  ion:,  over  successive  x  _  m  »  os  i(|  -  •;  in  , i( 

‘  i  pi  e  ement  bodies,  or  successive  source  or  ’  *>  I  »*  ’ 

l,l  >w’ng  v  lo,  it-,-  di'trihut  ion  ;  over  the  original 

1  ■  1  ■  ,  <  ont  inne s  utt t  i  I  a  g  i  •  n  d  i  f  f  e rence  c  r  it  er  ion  •  t  1  jv  - ’  1  *.  -  •  '  I  ' 

i  -  r*  t  .  1  xpet  ience  with  the  program  has  shown  r,, 

that  the  i  >  ‘mput  e,J  pressure  coefficients  for  t  hr 
.  -  *  v  d  and  third  iterations  usual  lv  agree  to 
within  0,0  1  over  n<>  ;t  -if  the  body.  Since  t  he 
result-.  usually  converge  in  an  osi  illatorv  manner. 
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where  —  *-  slnu, 
dn 


-  Upeos(iVi**lll) 


where  f  -  u/Ue  and  v/Ue  are  the  final  solutions 
oi  the  thin  boundary  layer  Equations  (1)  and  (2) 
us  ins  the  average  values  of  Ue  obtained  in  the 
last  two  iterations  of  the  viscous- invisc  id 
calculation.  The  invisc id  velocity  U  = 

Ux"  +  Ur*-)^  and  the  invisc  id  flow  angle  u*  * 

tan-1  (Ur/Ux)  are  computed  from  the  final 

potent ial -flow  computation  with  the  blowing  velocity 

distribution,  equation  (10),  on  the  original  body 

surface,  where  l>r  and  U  are  the  invisc  id  radial 

and  axial  velocities,  respectively. 


Ill*  A  Computation  Procedure  for  the  Parabolized 
Reynolds  Equations  in  Ax Isymmet r ic  Flow 
Us Ing  Streamline  Coordinates  and  the  k  -  r 
Turbulence  Model 


The  natural  coordinate  system  shown  in  Fig. 

1  Is  one  in  which  one  coordinate  s  lies  along 
the  streamline  and  the  other  two  n,9  are  normal 
to  the  streamline  (Liepmann  and  Roshko  {24}).  For 
ax  isymmet ric  flow,  =  0  where  .)/»)**  =  0  0  is  in 

the  azimuthal  direction  for  ax isymmet ric  flow.  In 
the  Reynolds-average  Nav ier-Stokes  equations,  the 
diffusion  terms  along  the  mean  streamline  direction 
are  usually  very  small  and  are  neglected,  the 
Reynolds  equations  with  such  an  approximation  are 


-  —  —  —  =  ic  and  di^  *  rUdn  ■  rU  (cos'j 
ds  R 

The  transformed  Equations  (14)  and  (IS)  i 
(x,  40  coordinate  system  are 


ip  *  1  ~  f,  r  ii'l 
ds  r  dn  L  6  <*nj 


It  3U  , 

l«  cos  t  —  *  - 

CO  Si 

■  -  H 

3x 

,x  L  J 

-  s  ***<.< 

il 

r2  v  U 

r 

9«  1 

e 

+  U  1- 

v  U  ~ 

3i p 

L 

e  3<p 

•  .2  lift 

IH  cosu  —  =  sin' 
3x 


ini  i-W-r  ui-  fP 
3x  [^,1  t)l i'  0 


-  ifi-(20) 

P  on 


The  turbulent  eddy  viscosity  vp  in  Equation 
(19)  is  modeled  by  the  two  k-L  equations 
governing  the  interaction  of  the  turbulence  with 
the  mean  flow.  The  turbulence  kinetic  energy  is 
defined  as  k  =  ^i7p/2,  and  the  t  urhul  enre 
energy  dissipation  is  defined  as  c=v( du i/ dxm) 
(<luj/^xm).  It  is  assumed  that  the  eddy  viscosity 
is  determined  by  dimensional  analysis  as 


The  k-k  equations  developed  by  Hanjalic  and 
launder  {251  can  he  written  in  the  (s,n)  coor¬ 
dinate  system  as 


2s  r  3r\  ^  J  ~  \  3n  ) 


where  i  .  the  tot.il  mean  velocity  along  the 
t  n  im!  im  direction  s,  n  is  the  distance  normal 
i  t  h»  mi  an  streamlines,  r  is  the  radial  distance 
rr.-n  t  h»-  b, njv  axis  x,  p  is  the  mean  pressure,  t 
i  •  tli*  angle  -it  the  streamline  with  the  x-axis, 

*  r,  is  the  molecular  kinematic  visro- 
it  ,  and  j  is  the  turbulent  eddy  viscosity.  A 
•  t  r<  am  !  niv  f  ion  .  j  .  defined  tor  ax  isymmet  r  ir 


(i) 


d  ,  r  !’  dn  (  n 

di.ewn  in  t  i  ,*  .  1,  the  tr.mstormaf  ion  of  the 
i  '  | '  i .  1 1  ion  .  ’ii  the  (.,  ,)  coordinate  system  can  be 
■'faired  with  the  aid  of  the  following 
r •  1  at  ion sh i ps 


.  •(  i  >()  >y  »0  ■,  ..  •( )  ‘O  .  1M. 

i  !  -  *•  r  I  -'.hi  i  (18) 

■n  •  v  in  cj  >n  *•;  *x 


All  the  turbulence  diffusion  terms  in  Equat  ions 
(22)  and  (21)  have  been  neglected  except  the  tern 
(u  '  -  v  *- )  *ti  /  *s_{  k/  1)  hi/  ds  in  the  turbulence 
production.  Hanjalic  and  Launder  2*>  recon- 
mend  the  retention  of  this  term  to  emphasize 
tin'  role  of  irrotat  ional  deformations  in  pro¬ 
moting  energy  transfer.  The  constants  P.  ,  i‘\, 

,  C) ,  y  and  o,  as  given  bv  Hanjal  ic  and 
launder  21)  i  are  0.09,  1.44,  1.90,  4.44,  1.0,  and 
1.10,  respectively.  The  t ran sf ormat ion  ot  the  k- 
<*r|(iat  i<>ns  (22)  ami  (21)  into  the  (x,c)  ccorit  in.d  c 
system  has  been  made  by  Hogan  22  using  equat  ions 
(lb),  (  17)  ,  and  (18) ,  o.g. 
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where  the  values  of  8,  vq>,  and  3u/9y  are  obtained 
from  the  final  iteration  of  the  simple  viscous- 
invlscid  interaction  computations.  The  lower 
boundary  in  the  computation  domain  is  the  stream 
surface  tj»  which  lies  along  the  body's  stir  face 
and  along  r«0  in  the  wake.  The  boundary  conditions 
on  <J>*0  in  the  wake  are  r*(i*0,  and  3U/3i|^»9p/3i^ 

3k/  *4,=  3i:/3k»0.  The  boundary  conditions  on  the  body 
('J'o  at  T^O)  are  r»rc,  U«0,  >i  -  hq  =  tan'l  Cdro/dx)  , 
and  the  condition  for  pressure  provided  by  equation 
(20) 


rll2  0»8  t  =  0  .  -  '£  .  sin  --E  -  rl- 

ax  3n  ax  >n 


Again,  the  assumption  of  »  =•  -  u  *v  ilU/Ijn  is  used 
to  obtain  the  boundary  conditions  for  •  and  k  as 
estimated  by  equations  (2b)  and  (27)  on  the  body 
('V  ~  0)  .  However,  the  values  of  t  and  k  are  not 
taken  at  the  wall  hut  at  a  small  distance  from 
the  wall  (usually  nu*/v  -  50)  where  the  inner 
region  of  the  mixing  length  and  eddy  viscosity 


i 


0.4  r 

o 


sin 

r 


and 


(28) 


[  fr 

I'li 


”  2 

r 

r“  v 

sin  i 

4  1'  -■ 

.>■  it 

<x 

)  aig 

«r 

■  III 

: 

? 

»  r  ,  - 

in  '  .»* 
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is  valid. 

An  iterative  numerical  marching  procedure  has 
been  developed  by  Hogan  ‘22:  to  solve  ecju.it  ions  (  1  ‘M  , 
(20),  (24),  and  (25)  within  the  partially  parabolic 
flow  assumption  that  the  inconsistent  pressure 
field  downstream  is  communicated  to  the  upstream. 

The  overall  numerical  procedure  of  Hogan  •22  is 
used  here;  however,  the  numerical  details  tar  soloing 
the  k—  equations  (24)  and  (25)  have  been  improved 
in  this  paper. 


iV.  Numerical  Results 


. '.pntution  ptmedure  starts  at  a  sla- 

*  hod%  where  the  boundary  layer  is  thin 

•  down st  ream  into  the  wake.  As  shown 

.  f  hi 1 1 in d, 1 1  conditions  for  the  values 

•  *  i  *  ju.it  ions  t  1 'O  and  (2<0  over  the 

•  *  •  •  j  >  and  along  a  cylindrical  stream 

.?  .  <  f  ' i* •  hound.if.  laver/wake  are 

! ,  ,  r  riot.  ■.  a  1  Me  .  oM  .i  ined  from  the 
o'  i  •  id  intera*  t  ion  eomputat  in. in 
c  r  ;  !!.  1  be  stream  surface  ,\  is 

1  :..r  if  lie*.  « nt  irelv  outside  the 

f  nr.  1  i\fi  and  wake  of  the  flow 

•  r  •  •  i<  ,  f  :.t  I  nt-  .  of  k  and  in 

•  •  and  '.“'l  ar**  /ere.  Furthermore, 

'  '  ’  <•  di  ..sip.it  ion  is  assumed  to 

’  r‘  .  «  pr  -dm  t  ion  inside  the  thin 

.  .  r  the  inlet  plane,  e.g.t  >  -  -u  *v  " 

<•  r<  .  the  distributions  ot  *  and  k 
>  *  1  r  h  -und ar  .  l.ivcr  i  an  he  est  imated  by 


Hsing  the  numerical  procedure  developed  in 
Sections  II  and  Til,  calculations  were  performed 
for  four  bodies  for  which  experimental  data  were 
available.  These  are  designated  as  DV'SKIM’  axi- 
svmcnet  r  i<*  Model  1  (Huang  et  al.  T  ),  HTNSklK* 

ax  i  svmmet  r  ic  Model  5  (Huang  et  al.  '«  ),  Model  A, 

and  Model  B  (l.yon  2  h ,  2  7  j )  .  Table  1  1  ists  the 

following  geometri'  and  flow  parameters  tor  each 
ho<!v:  the  length  of  the  hodv  I  ,  the  maximum 
radius  of  the  hodv  rmax,  the  upstream  flow  velo¬ 
city  I’  ,  and  tin'  hodv  Reynolds  number  k  .  No  flow 
o  e 

separation  on  t he  four  models  was  measured  or 

predicted.  The  measured  and  computed  mean  flow 
charaeterist  ics  over  the  sterns  ot  the  tour 
models  are  shown  in  Figures  2  through  h.  In 
all  eases,  the  wall  fund  ions  used  in  the  k- 
turhulence  model  are  taken  as  the  values  ot  k 
and  •  empaled  i rnm  equal  ions  (2h),  (27)  and 
(28)  with  nu^/y  set  equal  t o  50.  The  simple 
v  i  seou?.  -  inv  i  se  id  interact  ion  computat  ion  i  s 
the  romput  at  ion  procedure  out  l  ined  in  Sect  ion  I  1 
and  the  met  hod  designated  as  the  "para  ho  l  i.  ed 
N .  -S .  calcnlat  ion"  is  the  numerical  ptocedute 
of  solving  the  parabol  i/.ed  kevno  1  ds-avi -raged 
Nnv ier-Stokes  equations  using  streamline 
coordinates  and  the  k--  turbulence  model  sum- 
mar  iz ed  in  Seel  ion  Ill. 
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Norm.  1 1  Reynolds  Stress  Profiles  Across  Thick 
Stern  Houml.irv  J. ivers 


Figure  b  -  (Continued) 


- -  k~*  CALCULATIONS 

□  MEASUREMENTS  (3) 


k/U 


PARABOLIZED  N.-S.  CALCULATIONS  WITH  K-€ 


-  PARABOLIZED  N.-S.  CALCULATIONS  WITH  k-< 

- SIMPLE  VISCOUS- IN  VISCID  INTERACTION  COMPUTATIONS 

O  MEASUREMENTS  (LYON.  27 ( 


(b)  Distribution  of  Nondimensiona I 
Frictional  Velocity 


(c)  Distribution  of  Pressure  Coefficient 
on  the  Stern 


h/L  -  0  7000 


Figure  4  -  Comparison  of  the  Measured  and  Computed  Mean  Flow 
Characteristics  Over  the  Stern  of  Lyon  Model  A 


Figure  3  -  (Continued) 

-  PARABOLIZED  N.-S-  CALCULATIONS  WITH  k-t 

----  SIMPLE  VISCOUS- INVISCIO  INTERACTION  COMPUTATIONS 
O  MEASUREMENTS  (4) 


*/l  *  0  704  k/L  «  0  831  "  0.873 


-‘P 


»/L  -  0  987 


(e)  Distributions  of  Pressure  Coefficients 
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x/L  -  0.704 


x/L  -  0.831 


x/L  -  0.873 


NOTATION 


Cp  Coefficient  of  pressure  *  p  -  pQ 

1/2 pU02 

Cj  Numerical  constant  in  the  k-C  model  -  l.AA 

Cj  Numerical  constant  in  the  k-E.  model  -  1.90 

Cj  Numerical  constant  in  the  k-c  model  =  A.AA 

Cp~  C  Numerical  constant  is  the  k-C 

model  »  0.09 

k  Turbulent  kinetic  energy 
L  Length  of  the  body 

l  Mixing  length 

Mixing  length  in  the  inner  region  of 
boundary  layer  defined  in  equation  (A) 

Mixing  length  in  the  outer  region  of 
boundary  layer  defined  in  equation  (A) 

n  Distance  measured  normal  to  the  mean 
streamlines 

p  Mean  pressure 

pD  Upstream  pressure 

R  Radius  of  curvature  of  a  streamline 
Rc  Reynolds  number  of  the  flow  *  UQL/v 
r  Radial  distance  measured  from  body  axis 
rmax  ttexiam11'  radius  of  the  body 
rc  Body  radius 

s  Distance  measured  along  a  streamline  or 
body  meridian 

U  Total  mean  velocity  along  the  streamline 
d irect ion 

Inviscid  (edge)  velocity  used  for  thin 
boundary  layer  computation 

l!n  Free-stream  velocity 

Up  Total  inviscid  velocity  =  (U^x  + 

Uf  Inviscid  radial  velocity 

f'x  rnviscid  axial  velocity 

u  Mean  velocity  component  parallel  to  the 
hodv  meridian 

u  Turbulent  velocity  In  the  s-d Irect ion 

um  Modified  value  of  u  defined  in  equation 
(11) 

wx  Mean  axial  velocity 
u*  Friction  velocity  =  (rw/t-)5 


v  Mean  velocity  component  normal  to  the 
body  meridian 

v“*  Turbulent  velocity  in  the  y-  or  n- 
d irect ion 

vr  Mean  radial  velocity 

w  Turbulent  velocity  in  the  0  direction 

x  Coordinate  measured  along  the  axis  of 
the  body  measured  from  the  nose 

Beginning  x-station  of  the  partially 
parabolic  k-e,  calculations 

y  Distance  measured  from  the  body  surface 
normal  to  the  body  meridian 

yc  Value  of  y  at  which  -  v0  or  ■  l0 

a  Angle  of  streamline  with  the  x-axls 

otj  Angle  of  inviscid  streamline  with  the 
x-axis 

oto  Angle  of  body  surface  with  the  x-axis, 
a0  =  tan-1  (dr0/dx) 

^tr  Klebanoff's  interra it tency  factor  defined 
in  equation  (3) 

6  «  <$995  Boundary  layer  thickness  *  position 
where  the  velocity  is  99.5%  of  the 
potential  flow  velocity 

<Sr  Radial  distance  of  boundary  layer  thickness, 
6r  =  6cosa0. 

6*  Ax i symmetric  displacement  thickness  defined 
in  equat ion  (8) 

c  Turbulent  kinetic  energy  dissipation 
0  Azimuthal  angle 

k  longitudinal  curvature 

A  Displacement  area  defined  in  equation  (7) 

V  Molecular  kinematic  viscosity  of  the  fluid 

V  Total  effective  viscosity  =  V  +  Vf 

Inner  eddy  vicosity  defined  in  equation 
(3) 

vQ  Outer  eddy  viscosity  defined  in  equation 

O) 

\~y  Turbulent  eddv  viscosity 

♦  -  Constant  fluid  dens  it  v 

Turbulent  Prandt 1  number  used  in  the 
k  equat  ion  =  1.0 

\  Turbulent  Prandt I  number  used  in  the 
equat  ion  =  1.30 
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Extensive  measurements  of  the  pressure 
field  were  made  for  the  flows  past  Models  1  and 
5  (Huang  et  al.  [3,41.  Comparisons  of  the  com¬ 
puted  with  the  experimental  pressure  fields  for 
the  two  bodies  are  given  In  Figures  2e  through 
3e.  The  computed  pressure  profiles  and  the 
experimental  pressure  fields  for  Model  1  (Figure 
2e)  are  within  !• 5t;  for  Model  5  (Figure  3e )  the 
results  differ  by  2X.  The  computed  pressure 
profiles  agree  very  well  with  the  experimental 
results  above  the  boundary- layer  region. 

The  differences  between  the  computed  pressure 
coefficients  by  the  two  methods  are  all  less 
than  1%  outside  the  displacement  surface  and 
less  than  2Z  near  the  body.  This  indicates 
that  the  simple  viscous- inviscld  interaction 
computation  procedure  using  the  displacement 
body  concept  and  the  revised  mixing  length 


surface  and  is  about  two  percent  near  the  body. 
Only  near  the  very  tail  end  of  the  body  are  the 
parabolized  Nav for -Stokes  calculations  found  to 
model  the  flow  better  t  Ivin  t  lie  simple  method. 

The  prediet  ton  of  surface  pressure  coefficients 
and  fr'-tion  velocities,  normalized  axial  and 
radial  e. nettles  and  cross-stream  pressure  co¬ 
efficients  are  in  close  agreement  with  the  experi 
mental  data  for  four  models  having  attached  flows 
Except  tor  the  surface  friction  velocities  near 
the  tail  end  of  the  body,  an  overall  agreement  of 
t  lie  above  quantities  by  the  two  computation  method 
lias  been  obtained.  The  developed  simple  and  ef¬ 
ficient  v  i  scons-  inv  isc  id  interaction  computation 
procedure  ran  ho  used  as  a  design  tool  to  compute 
the  cross-stream  velocity  and  pressure  variat  ion 
across  the  thick  stern  boundary  layer  for  many 
practical  naval  applications. 


correctly  determined  the  essential  features  of 
the  thick  stern  boundary  layer.  Only  near  the 
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The  results  of  the  turbulent  kinetic  energy 
calculations  for  Model  1  and  Model  3,  together 
with  experimental  results,  are  given  in  Figure  6. 
The  agreement  of  the  computed  k  with  the  experi¬ 
mental  data  Is  good.  These  results  are  encouraging 
and  indicate  that  the  use  of  the  inner  mixing 
length  with  the  k-L  model  gives  a  good  approxima¬ 
tion  to  the  turbulent  field.  In  the  wake,  the 
agreement  is  satisfactory. 

The  measured  and  computed  Reynolds  stress 
profiles,  -u'v',  are  presented  for  Model  1  and 
Model  5  in  Figure  7.  The  agreement  of  the 
computed  results  and  the  experimental  data  is  also 
satisfactory  for  model  1.  The  measured  values  of 
-u  v'  for  model  5  are  higher  than  the  predicted 
values  at  X/L  >  0.87. 

Overall,  the  agreement  between  the  measured 
and  calculated  results  is  encouraging.  For  most 
of  the  ilow  field,  the  velocity,  and  the  pressure 
profiles  are  correctly  predicted  by  the  two 
methods.  The  wall  pressure  and  shear  stress 
distributions  computed  by  the  two  methods  also 
agree  well  with  the  experimental  data.  The 
measured  distributions  of  turbulence  kinetic 
rncrgv  k  and  Reynolds  stress  are  satisfactorily 
predicted  by  the  k-»  turbulence  model  used. 

V.  Conclusion 

Two  viscous- inv isc  id  interaction  computa¬ 
tion  procedures  are  presented.  One  method  solves 
the  parabolized  Reynolds-averaged  Nav  ier-Stokes 
equations  using  streamline  coordinates  and  t he 
k-  turbulence  model  and  t he  other  method  solves 
the  simpler  thin  boundary  layer  equations  using 
the  1 . ighthill  displacement  body  concept  and  the 
revised  mixing  length  for  the  thick  boundary  layer. 
I'he  tangent  ini  velocities  computed  bv  the  thin 
boundary  laver  equations  are  adjusted  to 
account  for  the  inv isc  id  influence  in  the 
simple  method.  The  maximum  difference  in  the 
computed  axial  and  radial  velocities  between  the 
two  methods  is  about  two  percent  of  t  fie  free¬ 
st  ream  ve  1  iw  it  v  and  the  maximum  difference  in  l  lit- 
computed  cross -st ream  pressure  coefficients  is 
lev.  than  one  percent  outside  the  displacement 
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TABLE  1  -  FLOW  AND  BODY  GEOMETRY  PARAMETERS 


L 

(m) 

rmax 

(m) 

(m/a) 

Re 

Model  1 

3.066 

0.1398 

30.48 

6.60  x  106 

Model  5 

2.910 

0.1398 

45.72 

9.30  x  106 

Model  A 

1.778 

0.1778 

17.88 

2.09  x  106 

Model  B 

1.778 

0.1778 

17.88 

2.05  x  106 

The  numerical  results  display  a  thickening  of  the 
turbulent  region  In  the  stern/wake  regions  of  the 
flows.  The  displacement  bodies  diverge  signif leant 2 y 
from  the  physical  bodies  near  the  stem  and  continue 
into  the  wake  with  slowly  decreasing  radii.  In 
Figures  2a  and  3a,  the  computed  displacement  body 
and  boundary  layer  thickness  are  compared  to  the 
values  of  <S*  and  <5  obtained  from  the  data  of 
Huang  et  al.  I  3,4  I.  For  both  bodies,  the 
computed  ft*  and  ft  1 ie  slightly  below  the  ex¬ 
perimental  results  in  the  stern/wake  region, 
hut,  overall,  the  agreement  with  the  experiments 
i  s  good . 

The  distributions  of  the  frictional  velo¬ 
city  u*  are  shown  in  Figures  2b,  3b,  4b  and  5b 
and  the  wall-pressure  coefficient  cp  are  shown 
in  Figures  2c,  3c,  4c,  and  5c  for  tne  four 
bodies.  The  pressure  distribution  computed  by 
the  parabolic  N.-S,  method  for  Model  1  (Figure  2c) 
has  a  large  trough  at  the  inflected  stern.  In  this 
region  of  the  body,  the  surface  and  the  stream¬ 
lines  near  the  surface  have  a  marked  change  In 
curvature.  As  the  streamlines  change  curvature 
from  convex  to  concave,  the  pressure  gradient 
changes  from  adverse  to  favorable.  Following  the 
concave  part  of  the  stern,  the  streamline  cur¬ 
vature  becomes  convex  again,  with  a  corresponding 
rise  in  the  pressure  on  the  wall.  The  computed 
wall  shear  stress,  given  by  pu^,  drops  rapidly 
in  the  adverse  pressure  gradient  region  of  the  flow. 
Accompanying  the  sharp  drop  in  the  wall  pressure, 
the  wall  shear  stress  rises  steeply.  With  the 
final  cliange  of  curvature  of  the  streamlines  at 
the  tail,  u*  drops  almost  to  zero.  The  computed 
pressure  distribution  agrees  well  with  the 
experimental  data,  liaving  a  maximum  percentage 

difference  of  1%  of  the  total  head  ill  /2. 

o 

However,  the  computed  u*  distrihnt  ion  bv  flu* 
par  a hoi  i/eti  method  ranches  .1  smaller  value  than  the 
experimental  data  near  the  tail  of  the  stern,  being 
.’V  ond>*r  the  e.xper  iroent  a  1  values  at  x/1  -  0.97. 

•  he  i  omput  ed  u*  and  rp  distributions  for  Model  4 
(figures  3h  anil  lv )  display  tile  same  t  vpe  of 
behavior  a.  exhibited  for  Model  1.  Fxcept  the  com¬ 
puted  ut  distribution  bv  the  simple  method  has  a 
higher  value  t ban  the  experimental  data  for  x/l 

For  Model  5,  the  agreement  with  the  experi¬ 
mental  data  is  good  for  both  the  wall  frictional 
velocity  and  the  wall  pressure  coef f  ic  ient .  As  is 
evident  bv  the  steep  drop  in  u*  in  the  adverse 
pressure  gradient  region  of  the  flow,  the  flow  about 
this  hodv  is  very  near  to  separation  at  x/l.  = 

'i.’M.  The  computed  u*  and  wall  cp  for  flow  past 
Model  A  are  given  in  Figures  3b  and  3c.  As 
opposed  to  Afterbodies  1  and  5,  the  stern  of  Model 
A  is  not  inflected.  Therefore,  over  the  aft  region 
of  the  body,  the  pressure  gradient  remains  adverse 
up  to  the  tail  of  the  hodv  and  u*  steadily  decreases 
to  zero.  The  computed  wall  pressure  distribution 


lies  slightly  above  the  experimental  pressure 
distributions  of  less  than  3£  at  x/l.  *  0.95.  The 
shape  of  Model  B  Is  character ized  by  a  sharply 
sloping  stern  region  (Figures  5a,  5b,  and  5c). 

The  values  of  the  surface  angles  near  90°at  the  stern 
1  caused  some  numerical  difficulty  for  the  parabolized 
(  N.-S.  computer  code.  However,  good  agreement 
of  the  computed  wall  pressure  with  the  experi¬ 
mental  data  is  obtained,  with  a  maximum  dif¬ 
ference  of  less  than  4%  occurring  at  x/l.  *  0.95 
(Figure  5c).  The  sharp  decrease  in  u* 

I  indicates  that  the  flow  is  nearing  separation 
i  as  the  body  sharply  turns  downward  at  the  tail 
(Figure  5b),  The  differences  between  the  com¬ 
puted  values  of  uA  and  cp  by  the  parabolized 
N.-S.  code  and  the  simple  viscous- invisc id  interac¬ 
tion  code  are  generally  small  except  near  the 
tall  end  of  the  stern. 

Figures  2d,  3d,  4d,  and  5d  present 
detailed  comparisons  of  the  computed  velocity 
fields  to  the  experimental  results  for  the  four 
bodies.  In  Figure  2d,  the  computed  velocity 
profiles  show  remarkable  agreement  with  the 
experimental  data  for  Model  1  up  to  the  tall  of 
the  body.  The  computed  velocities  near  the 
tail  and  in  the  wake  are  also  in  good  agreement 
with  the  experimental  results  (Figure  2d).  The 
largest  discrepancy  occurs  immediately  behind 
the  body  (x/I.  =  1.0076).  Farther  into  the  wake, 
the  agreement  with  experiment  is  very  good,  as 
is  evident  by  the  profiles  at  x/l.  *  1.1820. 

Computed  and  experimental  velocity  profiles  for 
Model  5  are  presented  In  Figure  3d.  The 
overall  agreement  with  the  experimental  data 
Is  good.  Figures  4d  and  5d  present  the  total 
velocity  profiles  for  the  flows  past  Models  A 
and  B.  The  agreement  of  the  computed  velocity 
fields  with  the  experimental  velocities  Is  very 
good  with  a  maximum  difference  of  less  than  27. 
Detailed  comparisons  of  the  parabolized  N.-S. 
calculations  and  the  simple  v iscous- inv isr id 
Interaction  computations  for  tlu*  four  bodies 
are  shown  in  Figures  2  through  5.  For  Model  l 
(Figures  2b  and  2r) ,  the  u*  and  cp  distributions 
are  in  general  agreement  up  to  x/l.  =  0.95.  The 
simple  method  does  not  predict  as  steep  a  drop 
in  u^  at  the  concave  part  of  the  hodv  as  do  the 
parabolized  N.-S.  calculations.  In  addition, 
there  Is  no  trough  at  the  stern  in  the  pressure 
distribution  computed  bv  the  simple  method  ( Figure 
2c).  The  velocity  profiles  computed  for  Model  l 
(Figure  2d)  from  the  two  methods  are  consistent  . 

For  Model  5,  the  parabolized  N.-S.  calculation 
correctly  predicts  the  steep  drop  in  u*  and  the 
pressure  trough  near  the  tail  end  of  the  stern 
(Figure  3b),  The  velocity  profile  comparisons 
for  Model  5  (Figure  Id)  demonstrate  again  the 
agreement  of  the  two  methods  (within  27),  The 
calculations  for  u*  and  rp  for  Model  A  (Figures 
4b  and  4c)  demonstrate  that  the  calculations  agree 
up  to  the  tail  of  the  hodv.  At  the  tail,  the  par¬ 
tially  parabolic  N.-S.  calculation  of  the  u*  dis¬ 
tribution  drops  rapidly  to  zero  as  opposed  to  the 
simple  method  which  shows  a  slight  upturn.  Both 
pressure  distributions  turn  downward  at  the  tail, 
with  the  parabolized  N.-S.  pressure  distribution 
reaching  a  slightly  higher  maximum.  The  computed 
velocity  profiles  for  Model  A  (Figure  4d)  are  again 
consistent.  The  comparisons  of  n*,  <p,  and  velocity 
for  Model  B  (Figures  4b  and  Sd)  show  the  same 
behavior  as  demonstrated  for  Model  A. 


c 


0  1* 


Figure  7  -  (Continued) 

-  fc-f  CALCULATIONS 

□  MEASUREMENTS  (4) 

k It  •  0.7031  1  1  >/L  *  0  ,7,7 


(b)  DTNSRDC  Axisymmetric  Model  5 


INITIAL  DISTRIBUTION 


Copies 

Copies 

1 

WES 

11  NAVSEA 

1 

SEA 

99612  (Library) 

1 

U.S.  ARMY  TRAS  R&D 

1 

SEA 

05R24  (J.  Sejd) 

Marine  Trans  Div 

1 

SEA 

55W3  (E.  Comstock) 

1 

SEA 

55W33  (W.  Sandburg) 

3 

ONR/432F  Whitehead, 

1 

SEA 

55W31  (W.  Louis) 

Lee,  Reischman 

1 

SEA 

55W31  (G.  Jones) 

1 

SEA 

55N2  (A.  Paladino) 

2 

NRL 

1 

SEA 

55N1  (S.G.  (Wieczorek) 

1  Code  2027 

1 

SEA 

63R31  (T.  Peirce) 

1  Code  2629 

1 

SEA 

56X12  (C.R.  Crockett) 

1 

SEA 

62R41  (L.  Pasiuk) 

1 

ONR/Boston 

1  NAVFAC/032C 

1  ONR/Chicago 

1  NADC 

1  ONR/New  York 


1 

NAVSHIPYD  PTSMH/Lib 

1 

ONR/Pasadena 

1 

NAVSHIPYD  PHILA/Lib 

1 

ONR/San  Francisco 

1 

NAVSHIPYD  NORVA/Lib 

1 

NORDA 

1 

NAVSHIPYD  CHASN/Lib 

3 

USNA 

l  Tech  Lib 

1 

NAVSHIPYD  LBEACH/Lib 

1  Nav  Sys  Eng  Dept 

1  B.  Johnson 

2 

NAVSHIPYD  MARE 

1  Lib 

3 

NAVPGSCOL 

1  Code  250 

1  Lib 

1  T.  Sarpkaya 

1 

NAVSHIPYD  PUGET/Lib 

1  J.  Miller 

1 

NAVSHIPYD  PEARL/Code  202.32 

1 

NOSC/Lib 

1 

NAVSEC,  NORVA/ 6660.03,  Blount 

1 

NCSC/712 

12 

DT1C 

1 

NCEL/131 

1 

AFOSR/NAM 

1 

NSWC,  White  Oak/Lib 

1 

AFFOL/FYS,  J.  Olsen 

1 

NSWC,  Dahlgren/Lib 

2 

MARAD 

1  NUSC/Lib  1  Div  of  Ship  R&D 

1  Lib 


21 


Copies 

1  NASA/HQ/Lib 

3  NASA/Ames  Res  Ctr,  Lib 

1  D.  Kwak 
1  J.L.  Steger 
1  Lib 

2  NASA/Langley  Res  Ctr 

1  Lib 

1  D.  Bushnell 
1  NBS/Lib 

1  NSF/Eng  Lib 

1  LC/Sci  &  Tech 

1  DOT/Lib  TAD-491.1 

2  MMA 

1  National  Maritime  Res  Ctr 
1  Lib 

3  U  of  Cal/Dept  Naval  Arch, 

Berkeley 
1  Lib 

1  W.  Webster 
1  R.  Yeung 

2  U  of  Cal.,  San  Diego 

1  A.T.  Ellis 
1  Scripps  Inst  Lib 

1  li  of  Cal.,  Santa  Barbara/Tulin 

4  CIT 

1  Aero  Lib 
1  T.Y.  Wu 
1  A.J.  Acosta 
1  D.  Coles 

1  California  State  University/ 

Cebeci 

1  Catholic  U  of  Amer/Civil  & 

Mech  Eng 

1  Colorado  State  U/Eng  Res  Ctr 

1  Cornell  U/Shen 


Copies 

2  Harvard  U 

1  G.  Carrier 
1  Gordon  McKay  Lib 

1  U  of  Illinois/J.  Robertson 

4  U  of  Iowa 

1  Lib 

1  L.  Landweber 

1  V.C.  Patel 
1  C.J.  Chen 

1  Johns  Hopkins  U/Lib 

4  MIT 

1  Lib 

1  J.R.  Kerwin 
1  T.F.  Ogilvie 
1  J.N.  Newman 

2  U  of  Minn/St.  Anthony  Falls 

1  Lib 
1  R.  Arndt 

1  U  of  Mich/NAME/Lib 

1  U  of  Notre  Dame/Eng  Lib 

1  New  York  U/Courant  Inst/Lib 

3  Penn  State 

1  B.R.  Parkin 

1  R.E.  Henderson 

1  ARL  Lib 

1  Princeton  IJ/Mellor 

1  U  of  Rhode  Island/F.M.  White 

1  Science  Application,  Inc. 

Annapolis,  MD 
C.  von  Kerczek 

1  SIT/Lib 

1  U  of  Texas/Arl  Lib 

1  Utah  State  U/Jeppson 


1 


Southwest  Res  Inst 
1  Applied  Mech  Rev 

Stanford  U 
1  Eng  Lib 

1  R.  Street,  Dept  Civil  Eng 
1  S.J.  Kline,  Dept  Mech  Eng 

Stanford  Res  Inst/Lib 

U  of  Virginia/Aero  Eng  Dept 

U  of  Washington/Arl  Tech  Lib 

U.S.  Naval  Academy 

VP  I 

1  Dept  Mech  Eng 
1  J.  Schetz,  Dept  Aero  & 
Ocean  Eng 

Webb  Inst 
1  Lib 
1  Ward 

Woods  Hole/Ocean  Eng 
Worchester  Pi/Tech  Lib 
SNAME/Tech  Lib 
Bell  Aerospace 

Bethlehem  Steel/Sparrows  Point 

National  Science  Foundation/ 
Eng  Div  Lib 

Bethlehem  Steel/New  York/Lib 

Boeing  Company/Seattle 
1  Marine  System 
1  P.  Rubbert 
1  G,  Paynter 
1  H.  Yoshihara 

Bolt,  Beranek  &  Newman/Lib 

Cambridge  Acoustical 
Associates,  Inc. 


Exxon,  NY/Design  Div/ 

Tank  Dept 

Exxon  Math  &  System,  Inc. 

1  Ceneral  Dynamics, 

EB/Boatwright 

1  Flow  Research 

1  Gibbs  &  Cox/Tech  Info 

2  Grumman  Aerospace  Corp 

1  Lib 

1  R.E.  Melnik 
1  Hydronautics/Lib 

1  Lockheed,  Sunnyvale/Waid 

1  Lockheed,  California/Lib 

1  Lockheed,  Georgia/Lib 

2  McDonnell  Douglas,  Long  Beach 

1  T.  Cebeci 
1  J.L.  Hess 

1  Newport  News  Shipbuilding/Lib 

1  Nielsen  Eng  &  Research 

1  Northrop  Corp/Aircraf t  Div 

1  Rand  Corp 

1  Rockwell  International 

1  B.  Ujihara 

1  Sperry  Rand/Tech  Lib 

1  Stanford  Research  Inst/Lib 

1  Sun  Shipbuilding/Chief 

Naval  Arch 

1  TRW  Systems  Group/Lib 

1  TRACOR 


Copies 

Copies 

Code 

Name 

1 

United  Technology/East 

1 

1843 

H.  Haussling 

Hartford,  Conn 

1 

19 

M.M.  Sevik 

2 

Westinghouse  Electric 

1 

1905.1 

W.  Blake 

1  M.S.  Macovsky 

1  Gulino 

1 

194 

F.S.  Archibald 

1 

Woods  Hole  Oceanographic,  Inc./ 

10 

5211.1 

Reports  Distribution 

Ocean  Eng  Dept 

1 

522.1 

TIC  (C) 

1 

522.2 

TIC  (A) 

CENTER  DISTRIBUTION 


Copies 

Code 

Name 

1 

1500 

W.B.  Morgan 

1 

1504 

V.J.  Monacella 

1 

1506 

S .  Hawkins 

1 

1508 

R.  Boswell 

1 

152 

W.C.  Lin 

1 

1521 

W.  Day 

o  1 

1522 

C.C.  Hsu 

1 

1522 

M.B.  Wilson 

1 

1522 

C.H.  Sung 

1 

154 

J.  McCarthy 

1 

154 

P.  Granville 

1 

1540.1 

B.  Yim 

1 

1540.2 

R.  Cumming 

30 

1542 

T . T .  Huang 

1 

1542 

N.C.  Groves 

1 

1542 

G.S.  Belt 

1 

1542 

Y.T.  Lee 

1 

1542 

M.S.  Chang 

1 

1542 

R.W.  Burke 

1 

1544 

F.  Peterson 

1 

156 

D.S.  Cieslowski 

1 

1561 

G.  Cox 

1 

1564 

J.  Feldman 

1 

1606 

T.C.  Tai 

1 

1802.1 

H.  Lugt 

1 

1840 

J.  Schot 

24 


END 

FILMED 

5-85 


DTIC 


